Charge Fluctuations from the Chiral Magnetic Effect in Nuclear Collisions 
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We derive a nonlocal effective Lagrangian for the chiral magnetic effect. An electric field is 
generated by winding number fluctuations of the nonabelian gauge field in the presence of a strong 
magnetic field. We estimate the magnitude of charge asymmetry fluctuations with respect to the 
reaction plane induced by the chiral magnetic effect in relativistic heavy ion collisions to be less 
than 10 -6 , several orders of magnitude smaller than the signal observed in the STAR experiment. 
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As pointed out by Fukushima, Kharzeev, and War- 
ringa [l|, a magnetic field induces a parallel electric cur- 
rent in the presence of a winding number transition in 
the color-SU(3) gauge field. This phenomenon, which 
has been called the chiral magnetic effect, can in princi- 
ple be studied in off-central relativistic heavy ion colli- 
sions. Since the magnetic field is oriented parallel to the 
reaction plane, fluctuations of the winding number of the 
gauge field are expected to cause fluctuations in the elec- 
tric charge distribution of emitted hadrons with respect 
to the reaction plane. Possible indications for such a 
phenomenon have been found in Au+Au collisions at the 
Relativistic Heavy Ion Collider, but the origin of the ob- 
served charge fluctuations among the final-state hadrons 
is still under debate [Jj. A more reliable estimate of the 
magnitude of the fluctuations caused by the chiral mag- 
netic effect would therefore be of great interest. Here we 
attempt such an estimate. 

The electric current induced by an external magnetic 
field in the presence of a winding number transition can 
be expressed in the form 



cr x B, 



(1) 



where a x is known as the chiral magnetic conductivity 
Q. Kharzeev et al. gave the following expression for a x 
in a thermal quark-gluon plasma [l[ : 



E 



16tt 4 T 2 



J alt (E a -B a ), 



(2) 



where E a ,B a denote the SU(3) gauge fields, g is the 
SU(3) gauge coupling, and e/ is the electric charge of 
each quark flavor. We proceed by redcriving in a 
slightly more general context and construct an effective 
infrared Lagrangian encoding the chiral magnetic effect. 

It is useful to first remind ourselves of the derivation 
of the electric analogue of eq. (jSJ, i. e. Ohm's law j = crE 
by starting from the general expression for the current 
density: 



(3) 



where n is the density of (nonrelativistic) charged parti- 
cles. An external electric field accelerates the particles, 



and they acquire the momentum 

Ap(i) = mv(t) = e [ dt'F,(t'). 
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(4) 



Denoting the average time between collisions that ran- 
domize the particle momentum as t c , the average addi- 
tional collective particle velocity is 
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We thus obtain for the collective current induced by the 
electric field: 
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E = o-E. 



(6) 



which is the Drude formula for the electric conductivity. 
In order to derive an effective action from this relation, 
we use the gauge A = — £E, dA = —tdE to obtain: 



C eS = - J(j)-dA = t J aE-dE=|^ dt'-E 2 



(7) 



The nonlocal effective Lagrangian becomes more familiar 
when expressed in terms of the Fourier transformed field 
E(cj): 



£ off = £|EH| 2 , 



(8) 



which shows that a /to is the imaginary part of the dielec- 
tric function e(uS). 

We now return to the chiral magnetic current ([T]). In 
order to derive an expression for a x , we need to split 
the particle content into left-handed and right-handed 
components: 
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The change in the difference between left-and right- 
handed particle densities is governed by the axial 
anomaly: 
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In order to calculate the dependence of (j) on the mag- 
netic field B and the excess of left-handed particles, we 
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introduce a chemical potential (i§ for the chirality. The 
energy of a quark whose spin is aligned (anti-aligned) 
with the magnetic field in the z-direction is 



(17) 



E(p,s z = +1/2) = ^ + to 2 + (2n + 2s z + l)eB 

= \/p 2 + m 2 + 2eB, 

E(p,s z = -1/2) = vV+w 2 . (11) 

where we have substituted the term associated to the 
individual Landau orbits (characterized by n) by p^ 2 and 
neglected the anomalous magnetic moment contribution. 
We will see that only the term linear in 2eB contributes, 
which can be calculated assuming that it is infinitcsimally 
small. This allows to substitute the discrete values of 
by continuous ones. The chiral chemical potential 
/i5 couples to chirality a. For a Landau level moving 
with momentum p z , the chirality is given by by a = 
2s z sgn(p 3 ). Note that while p_L changes with time, the 
helicity of a massless quark remains unchanged due to 
the vector nature of the electromagnetic field. 

We study the massless case in which the chiral charge 
density and the chiral magnetic current can be expressed 
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E[p,asgn(p 3 )} - afx 5 

x(l-2n F (E\p,asgn(p 3 )})) (12) 
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E\p,asgn(p 3 )] +<7(i 5 
x(l-2n F (E\p,asgn(p s )})), (13) 

where n F (E) = (e E / T + 1) _1 is the Fermi distribution. 

The leading divergences in these expression all cancel 
and one obtains simple finite results. The evaluation of 
the integrals is discussed in detail for the more difficult 
case (j) in Appendix A. The results are 

M5T 2 
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We can use these two equations to eliminate fi§ and to 
obtain the desired relation between the chiral quark num- 
ber density and the electromagnetic current: 

3e 2 B 
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J dt (E a ■ B a ). 



Different mechanisms provide three different upper 
limits to the time integrals in (|17|) : 

1. The transition time r sp h between gauge field con- 
figurations with different winding number; 

2. the time t x between many-body interactions that 
change the helicity of a quark; and 

3. the life-time tb of the strong magnetic field. 

In a quark-gluon plasma or a similar highly excited QCD 
medium without spontaneous chiral symmetry breaking, 
the transition time between states of different winding 
number is determine by the rate of so-called sphaleron 
transitions. These are diffusion processes of the gauge 
field across the semiclassical saddle point between con- 
figurations with different winding number. The char- 
acteristic time scale for these transitions is r sp h ~ 
[a^THl/as)}- 1 @. 

Helicity changing processes that are not related to 
sphaleron transitions are either governed by the nonzero 
current quark mass or by instantons, i. e. quantum fluc- 
tuations of the winding number of the gauge field. Per- 
turbativc interactions among massless quarks and be- 
tween quarks and gluons are helicity conserving. Cur- 
rent quark mass effects are negligible for u and d quarks 
and occur on a time scale of 1 fm/c for strange quarks. 
Instantons, i. e. topological tunneling processes, which 
occur on a time scale of 0.3 fm/c in the vacuum, are 
strongly suppressed at nonzero temperature. The helic- 
ity decay time r x therefore is expected to be much larger 
than the other time scales, which are estimated to be 
below 1 fm/c. 

The life-time of the strong magnetic field in a relativis- 
tic heavy ion collision at central rapidity is controlled by 
the longitudinal extent of the Lorentz contracted nuclei, 
which is tb ~ R/l, where R is the nuclear radius and 
(14) 7 = E cm /M]y is the Lorentz factor. Here E cm is the beam 
energy per nucleon in the center-of-mass frame and M/v is 
the nucleon mass. At the highest RHIC energy, 7 « 100; 
at the Large Hadron Collider (LHC) the factor will be 
7 w 3,500. 

Even at RHIC energies, tb will be comparable, if not 
shorter than the sphaleron transition time r sp h. This has 
two implications. One is that, in effect, the chiral mag- 
netic effect will be mainly sensitive to the time averaged 
magnetic field pulse 



(15) 



„ B / dt(E a -B a ). (16) 
16tt 4 T 2 

which agrees with the result obtained in ref. [l[ . Follow- 
ing the same argument as in the case of Ohm's law, we 
can now cast our result in terms of a nonlocal effective 
Lagrangian for the chiral magnetic effect: 
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dtB(t'). 



(18) 



The other is that the estimate of the chiral magnetic 
effective action, which assumed an approximately con- 
stant magnetic field, becomes questionable. To conclude 
this section, we write the effective action for the chiral 
magnetic effect, including the contributions from all light 
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quark flavors, as 
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(E-B) f dt" (E Q -B a ). (19) 
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A local electric current induced by the magnetic field 
during the early stage of a relativistic heavy ion collision 
only becomes observable when it is reflected in the final 
momentum distribution of charged particles. Asakawa 
et al. showed that the conversion of the local cur- 
rent into a final-state observable can occur owing to the 
collective outward flow of the matter. They related the 
charge asymmetry fluctuations with respect to the reac- 
tion plane to the vertical gradient of the expansion ve- 
locity at the time of particle freeze-out. Here we are 
pursuing a slightly different approach by considering the 
back reaction of the electric field to the chiral magnetic 
effect. The effective action (flU)) is linear in the electric 
field strength E. This means that the lowest total energy 
density is attained for a nonzero value of the electric field. 
To see this, we add f(T9|) to the free electric field energy 
density in the presence of the magnetic field: 



£he(E) = -/(B)E 2 . 



(20) 



where the function /(B) includes nonlinear quantum 
electrodynamic (QED) effects @. By explicit calculation 
one finds that for a constant magnetic field of magnitude 
B = (lOOMeV 2 ) these effects give rise to a very large 
/(B) sa 10, i. e. nonlinear QED effects can reduce the 
induced electric field strength substantially. Therefore, 
for any quantitative study one would have to calculate 
the strong QED effects for a realistic, time- and space- 
dependent magnetic field, which is a very demanding task 
indeed. The minimum of the field energy is obtained for 

e 2 o 2 rn r TB 
Emin = , { , x • B / dt" (E Q • B Q ) (21) 
8tt 4 T 2 /(B) J V ) \ ) 

We can also understand this relation as consequence of 
Maxwell's equation d~E/dt = — j. The nonvanishing elec- 
tric field accelerates quarks of opposite charge in opposite 
directions. As a result, the momentum space distribu- 
tions of positive and negative charged particles will be 
slightly shifted with respect to each other by the mo- 
mentum 

ApwgE min r B . (22) 

This momentum shift further deforms the thermal spec- 
trum of the emitted particles, which is already influenced 
by the collective flow v. The resulting spectrum 

d 3 N 

E p — aexp(- p n [Ep - v • (p Ap)]) , (23) 

where (3f = TJ 1 is the inverse freeze-out temperature 
and 7 = (1 — tj 2 ) -1 / 2 is the Lorentz factor associated with 
the collective flow, implies a charge ratio (with q = ±e): 

« exp (2er B v • E min /7» . (24) 



The sign of v • E m ; n fluctuates randomly between different 
coherence domains of the Chern-Simons number density 
E a • B a . The total charge asymmetry is thus given by 
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(25) 



where N denotes the number of independent domains. 

This formula assumes that the momentum kick re- 
ceived by the quarks during the action of the chiral mag- 
netic effect survives until freeze-out. One would expect 
that diffusive processes in the quark-gluon plasma result 
in a reduction of the initially generated momentum dif- 
ference between quarks with different charges. On the 
other hand, the small initial spatial separation between 
positive and negatively charged quarks may amplify the 
final momentum shift owing to the gradient of the col- 
lective velocity field in the directionperpcndicular to the 
reaction plane, as discussed in ref. [7|. 

We now estimate the magnitude of the expected fluc- 
tuating charge asymmetry. Using the result for the time 
integrated magnetic field in Au+Au collisions given in 
the Appendix of ref. Q , we have 
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where R, and Z are the nuclear radius and charge, re- 
spectively, and b is the impact parameter. Summing over 
up, down, and strange quarks and neglecting the quark 
masses, we find 
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(27) 



where |Qs|/V denotes the volume density of winding 
number transitions. Lacking a more reliable estimate, 
we assume one winding number transition per volume 
p 3 with p sa 0.5 fm during the presence of the strong 
magnetic field. Finally, we choose T pa 350 MeV as a 
proxy for the quark density during the strongly magne- 
tized phase of the heavy ion collision, and set /(B) sa 10. 
For Z = 79 and R = 7 fm this yields 



|eE n 



(6 MeV) 2 (b/R). 



(28) 



Finally, using the estimate tb ~ 2R/ r y cm with 7 cm = 100 
for the duration of the strong magnetic field, a final col- 
lective flow velocity |v| sa 0.5c, a freeze-out temperature 
Tf sa 150 MeV, and N = nR 2 /p 2 600 independent 
winding number domains, we obtain 



3.5 x 10 



_ 7 b_ 
R' 



(29) 



which is much smaller than the charge correlation ob- 
served in the STAR experiment 0, E| . 

This estimate can be improved in many ways. The 
most important improvement would be a more reliable es- 
timate of the winding number fluctuation density during 
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the strongly magnetized phase of the heavy ion reaction. 
This requires calculations of the contribution to Q5/V 
from the interaction between the color glass condensate 
clouds contained in the two colliding nuclei Q and of the 
sphaleron transition rate in pre-equilibrated quark-gluon 
plasma Q. Another approach to this problem involve 
lattice QCD simulations of the winding number fluctua- 
tion rate in the presence of a strong external magnetic 
field QMi- 

A second possible improvement concern the calcula- 
tion of the electric current induced by a winding number 
transition in the magnetic field, since the quark distribu- 
tion is unlikely to be thermalized during the early period 
when the magnetic field is strong. One way in which this 
estimate can be improved is by calculating the quark- 
antiquark distribution liberated by the interacting color 
glass condensates fl3j . 

Another important improvement would be the inclu- 
sion of charge transport processes during the hydro- 
dynamic expansion of the quark-gluon plasma and the 
hadronic freeze-out phase. Such calculations would facil- 
itate more realistic estimates of the final charge asymme- 
try fluctuations with respect to the reaction plane. 

Finally, the nonlinear QED effects have to be calcu- 
lated for a more realistic field configuration. Although 
is is not relevant for the scenario of "small" winding 
number coherence domains considered here, we mention 
a QED effect that may limit the chiral magnetic effect 
for large coherence domains. Electron-positron pair pro- 
duction limits the strength of stable electric fields to 
eE < eE max = 2m e /p, where m e is the electron mass 
and p denotes the winding number domain size. 

In summary, we have derived a nonlocal effective La- 
grangian for the chiral magnetic effect. The charge cur- 
rent corresponding to this Lagrangian agrees with earlier 
calculations. We show that a fluctuating electric field 
is generated by winding number fluctuations of the non- 
abelian gauge field in the presence of a strong magnetic 
field, and we estimated its effect on the charge asymme- 
try of the spectra of emitted hadrons. We estimated the 
fluctuating relative charge particle asymmetry in the fi- 
nal state and identified several opportunities for future 
improvements of the estimate. 
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Appendix A: Chiral electric current 

In order to calculate the dependence of (j) on the mag- 
netic field B and on the excess of right-handed parti- 
cles, we follow ref. [l[ and introduce a chemical potential 
/i,5 <C T for the chirality a = ±1. The energy of a quark 
whose spin is aligned (anti-aligned) with the magnetic 
field is 



E{p,s z ) = ^JpI + m 2 + (2n + 2s z + l)eB 
= v/p 2 + m 2 + (s + l)eB 
= E PtS . (Al) 

where s = 2s z = ±1 is twice the eigenvalue of the spin 
component along the magnetic field, here taken as the 
z-dircction. The eigenvalues of the transverse kinetic en- 
ergy are given by the Landau levels (pj_) n = 2neB. We 
will assume in the following that the magnetic fields is 
weak (eB <C T 2 ) so that the Landau orbits form a quasi- 
continuum and the sum over n can be approximated by 
an integral over pj_. 

The interaction of the quark with the magnetic field 
does not violate chirality. This implies the relation 
a = ssgn(p z ). We thus obtain the following distribu- 
tion functions for right- and left-handed quarks in the 
presence of a magnetic field: 

n a (p) = n F (E p . s - (j, 5 a), (A2) 

where n-p(E) denotes the Fermi-Dirac distribution. We 
thus find the following expressions for the net chirality 
and the average velocity density in the presence of the 
magnetic field: 

U = _ e y [ *P JbL. tanh ^ - . (A3) 
^ J (2tt) 3 E p , s 2T 1 ; 

We now expand for small p,§ and B up to linear order in 
[15 and the magnetic field B: 



d 3 p dE p ^ s d 



= -p 5 e(eB) 
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where the sum over chiralities has contributed a factor 2, This agrees with the result obtained by Fukushima et al. 

and Ep = p 2 + m 2 . In the following we will set m = 0. 

As (s + l)eB contributes only for s = +1 = crsgn(p z ), the 

integral over p z only runs over positive (negative) p z for 

positive (negative) chiralities a. Performing the angular 

integrals and setting E p = p = |p| for massless quarks, 

we obtain 

, z = -^/W#(^tanhA 
J 8tt 2 J yy dp \pdp 2T 



87r 2 J ' dx Vxcosh 2 
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